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Using first-principles plane- wave calculations study of electronic and magnetic properties of hypo- 
thetical two-dimensional structure of Li2N compound have been conducted. Calculations show, that 
electronic properties of this this structure can be inflenced by hydrogenation, which may change the 
system from wide-gap semiconductor to metal. Also, non-zero magnetic moment, equal to 1 can 
be generated by intruduction of H vacanies in hydrogenated structure. 

PACS numbers: 73.22.-f, 75.75.-c 



I. INTRODUCTION 

Since its discovery in 2004 graphene [l[ draws 
much attention because of unique features of this two- 
dimensional system. Graphene is composed of a sp^- 
bondcd carbon atoms forming honeycomb structure. It 
became famous for its very interesting electronic struc- 
ture with characteristic, linear energy dispersion near K 
point of Brillouin zone and many other features Q . 
Shortly after, experimental techniques allowed fabrica- 
tion of other new two-dimensional materials, like BN and 
M0S2 honeycomb structures 0] . The discovery of such 
stable two-dimensional materials triggered search for sim- 
ilar structures made from different compounds. Up to 
now many of these hypothetical structures constructed 
from silanene (2D Si) and germanene (2D Ge) 0, i, 
III-V compounds |H, SiC or ZnO i have been 
studied theoretically. Also, calculations show that 
graphene-like type of structure is not the only one possi- 
ble for two-dimensional material. This new class of boron 
sheets, composed of triangular and hexagonal motifs can 
be stabilized by interplay of three- and two-center bond- 
ing scheme [lO|. Another example of triangular sheet 
could be found in already known material, which is LiaN 
in its a phase. 

LiaN is a bulk material known to be a fast ion conduc- 
tor [s^]- LiaN is also known as a candidate for hydro- 
gen storage material due its high theoretical H2 capacity 
|13| . Bulk LiaN crystallizes in hexagonal structure which 
is characterized by P6 / mmm symmetry group, each ni- 
trogen atom is surrounded by eight lithium atoms. It has 
layered structure, one layer is Li2N and the other is of 
Li atoms only. Previous theoretical studies confirm ionic 
nature of bonding in this compound [ll|, [l^] ■ Since N- 
containing layer is rather weakly bound with two Li-only 
layers, it would be interesting to study electronic proper- 
ties of such two-dimensional structure (2DLi2N) - Fig[T^. 
Since this structure would have N atoms with dangling 
bonds, it would give opportunity to study influence of 
different atoms addition on them. For example addition 
of hydrogen atoms in case of graphene resulted in new 
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material which is graphane jl5| . 

Graphene and other nano-scale materials are recognized 
as future building blocks of new electronics technologies 
[l6j , including spintronics [13] • In the case of low (one- 
and two-) dimensional structures problem arises because 
of famous Mcrmin- Wagner theorem [Tsj , which prevents 
ferro- or antifcrromagnetic order to occur in finite tem- 
peratures, which is essential for practical application. 
This started the theoretical and experimental search for 
magnetism in graphene and other two-dimensional struc- 
tures. One of the most promising directions is emergence 
of magnetism in such structures as an effect of presence 
of local defects [l^ . According to works of Palacios et al. 
[20j and, independently, of Yazyev [2l| single-atom de- 
fects can induce ferromagnetism in graphene based mate- 
rials. In both cases, the magnetic order arises as an effect 
of presence of single-atom defects in combination with 
a sublattice discriminating mechanism. In the case of 
2DLi2Nrole of such defect could play non- hydrogenated 
N atom in hydrogenated structure. It would be then 
instructive to check influence of hydrogenation level on 
magnetic moment of the structure. 

In this paper electronic and magnetic structure of pure 
and hydrogenated 2DLi2N have been analyzed by means 
of ab-initio calculations. 



II. COMPUTATIONAL DETAILS 

To investigate electronic and magnetic properties of 
two-dimensional LiaN structures a series of ab-initio 
calculations have been conducted with use of DFT 
VASP code m, with PAW potentials For 
both spin-unpolarized and spin-polarized cases exchange- 
correlation potential has been approximated by general- 
ized gradient approximation (GGA) using PW91 func- 
tional [15. Kinetic energy cutoff of 500 cV for plane- 
wave basis set has been used. In all cases for self- 
consistent structure optimizations, the Brillouin zone 
(BZ) was sampled by 20 x 20 x 1 special k points. 
All structures have been optimized for both, spin- 
unpolarized and spin-polarized cases unless Feynman- 
Hellman forces acting on each atom become smaller than 
10-4 eV/A. A vacuum spacmg of 12 A was applied to 
hinder the interactions between 2DLi2Nmonolayers in ad- 
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jaccnt cells, (dop. kicdy supercell i jak liczone magn.) 
Bandstructure and density of states (DoS) calculations 
have been confirmed by use of WIEN2k code [2^ which 
implements the full-potential linearized augmented plane 
wave (FLAPW) method [13] • In this case for exchange 
and correlation generalized gradient approximation was 
used in the Perdew-Burke-Ernzerhoff (PBE) parameter- 
ization f3l|. 



III. RESULTS - ELECTRONIC STRUCTURE 

To study electronic properties of 2DLi2N, at first com- 
parison has been made with bulk material. For both cases 
lattice constants have been determined by total energy 
calculations and are found to be equal to 3.65 A for bulk 
(experimental value 3.63 A) and 3.57 A for 2DLi2N. In 
agreement with bulk LiaN is a semiconductor with 
non-direct bandgap equal to 1.15 eV between A (valence 
band) and F (conduction band) points. In contradiction 
to this, 2DLi2N has metallic nature. 
Two-dimensional structure is rather weakly bound - 
binding energy (defined as Ei, = Eat - Esheet where Eat 
is the energy of isolated atom(s) and Esheet is the total 
energy of two-dimensional structure) is equal to 10.36 eV, 
while binding energy of bulk structure is equal to 14.25 
eV. Also, two dimensional sheet would have N atoms with 
dangling bonds, such structure would be then rather un- 
stable with respect to foreign atoms addition. Graphane 
case suggests that it would be instructive to examine in- 
fluenece of hydrogenation on electronic structure in such 
cases as well as addition of lithium atoms. 
The nature of Li-N bond is ionic, as it can be seen from 
Fig. [5] showing charge density projected on [110] plane. 
Since every bond has both ionic and covalent character 
the level of ionicity can be estimated using difference 
between electronegativities of bonded atoms HI]. In 
the case of Li-N bond this difference equal to 2 suggests, 
that the bond is about 65% ionic and 35%covalent. This 
fact together with rather large lattice constant suggest 
that the structure of two-dimensional 2DLi2N can be 
low-bucklcd (LB) rather than plane (PL), according to 
puckering mechanism described in . To check this the 
series of calculations has been done, each with different 
distance in z direction between Li atoms and the plane 
on which N atoms lie (Az). The structure with minimal 
energy has been then optimized. Calculations show, that 
the buckled structure with Az = X lies 0.54 cV lower 
that the plane, which means that indeed the puckering 
mechanism stabilizes the structure. Both, plane and low- 
buckled structures can be seen on Fig [TJ 
Four structures have been then studied in two confor- 
mations, plane and low-buckled — two (PL and LB) 
with single H atom attached on top of each N atom 
(2DLi2N-|-H), two with two H atoms attached on both 
sides of N (2DLi2N-|-2H), two with single Li atom at- 
tached on top of each N atom (2DLi2N-|-Li), and two with 
two Li atoms attached on both sides of N (2DLi2N+2Li). 



Calculated binding energies, hydrogen-addition ener- 
gies, lithium-addition energies, and structure parameters 
are shown in Tab. 1 (plane structure) and Tab. 2 (low- 
buckled structure). In Tab. 1 d-N-H is the distance be- 
tween N and passivated H atoms, d-N-Li is the distance 
between N and passivated Li atoms, and Az is distance in 
the z direction between Li atoms and the plane on which 
N atoms lie. In tab. 2 Et is total energy of the structure 
with reference to plane structure. As can be seen, in all 
cases low-buckled structure is lower in energy, from 0.01 
eV for 2DLi2N-H2Li to 1.08 eV for 2DLi2N-h2H. So, from 
now on, all properties are referred to LB structure, unless 
stated otherwise. 

Binding energy of 2DLi2N+H, equal to 16.27 
eV, is comparable to binding energy of bulk struc- 
ture (both have the same number of atoms in the 
unit cell), which suggests that H-passivated two- 
dimensional structure would be no less stable than 
the bulk. Addition energies are defined as Eadd = 
i?tot (layer -faddatom) — iJtot (layer) — i?tot (addatom). 
High value of addition energy for single H atom equal to 
5.55 eV suggests, that hydrogen addition may stabilize 
the structure, making it more bound. Addition of second 
H atom to already H-passivated structure requires only 
1.44 eV, which means that 2DLi2N-|-2H structure is 
less stable than 2DLi2N-|-H. Addition energy of 2.11 
eV for single Li atom on pure layer is almost equal 
to addition energy for single Li atom on Li-passivated 
structure (1.74 eV), but such structures are less bound 
that H-passivated cases. 

MW -> 

The electronic structure of a single-layer 2DLi2N (not 
taking into account the Fermi level) exhibits typical 
structure of a semiconductor, with valence and conduc- 
tion bands separated by a band gap of 4 eV. Its metal- 
lic character is purely due to the position of Fermi level 
which is lowered by -2 eV from the middle of the band 
gap to the upper part of the valence band. This may be 
explained by a large contribution to the density of states 
of the valence band coming from the p-electrons from 
N atoms. The conduction band of a single-layer 2DLi2N 
originates mainly from the density of states of p-electrons 
of lithium. 

The relative positions of bands in band structure 
of 2DLi2N+H with single H atoms added, remain un- 
changed with respect to 2DLi2N-|-H with the exception 
of a single nitrogen p-electrons band, which is lowered by 
2 eV. The Fermi level is rised by 2 eV from the conduc- 
tion band to the middle of the band gap, and therefore 
2DLi2N-|-H is semiconductor. The rising of the Fermi 
level is caused by an increased contribution to the den- 
sity of states from lithium p-electrons (see Fig. 2e). The 
contributions of s-states originating from additional hy- 
drogen is insignificant. 

The same effect can be observed when the second H 
atom is added to the 2DLi2N-|-H layer. The Fermi level 
is rised again by 2 eV with respect to the position of 
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bands, and is placed at the bottom of the conduction 
band. Therefore 2DLi2N+2H is again metaUic. 

We can summarize the effect of H atom addition: each 
H atoms rises the Fermi level by 2 eV and the band struc- 
ture remains unchanged. 

Addition of lithium atoms to the single-layer 2DLi2N 
leads to quite different effects. The band structure of 
2DLi2N-|-Li is purely metallic with no band gap. The en- 
ergy band of p-clcctrons for additional litium atom con- 
nects the conduction band of the original 2DLi2N band 
structure at K and M points of the Brillouine zone with 
the valence band at T point, effectively nuUyfing the en- 
ergy gap. There exists non-zero density of states of p- 
electrons for additional litium atom directly at the Fermi 
level. The Fermi level is practicaly unchanged with re- 
spect to the band structure and remains at the top of the 
conduction band of the original 2DLi2N band structure. 

Addition of second lithium atom repeats the effects of 
the first one. There appears a new steep band of lithium 
p-electrons and the resulting structure remains metallic. 
The Fermi level is shifted up by about 1 eV with respect 
to the band structure, however is does not influence the 
metallic character of the compound. 

We can summarize the effect of Li atom additions: 
each Li atoms adds an energy band directly across the 
energy gap of the the original 2DLi2N band structure. 

MW 4- 



IV. CONCLUSIONS 

Ab-initio calculations have been conducted for hy- 
pothetical two-dimensional material 2DLi2N to investi- 
gate electronic and magnetic properties. Calculations 
show, that structure is much more stable when dangling 
bonds of nitrogen atoms arc functionalizcd with hydro- 
gen atoms. This hydrogcnation has very strong influence 
on on bandstructure, changing it from wide-gap semicon- 
ductor to metal. 

Magnetic properties are also interesting. In analogy to 
graphcnc and other two-dimensional materials it is possi- 
ble to generate non-zero magnetic moment by introduc- 
tion of distorsion. In the case of 2DLi2N the distorsion 
would be a two-hydrogen or hydrohen-lithium vacancy 
around the same nitrogen atom. This generates magnetic 
moment of I ■ Since bulk Li3N material has ususally 1- 
2% Li vacancies in Li2N layers [lj| such two-dimensional 
hydrogenated sheet would be almost naturally magnetic. 
These results may give a hint for experimentalists seek- 
ing for two-dimensional (magnetic) materials, which 
would be interesitng addition to growing family of two- 
dimensional materials. 
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FIG. 1. Atomic structure of plane (a) vs. low-buckled (b) 
2DLi2N (color online) 




FIG. 2. Valence charge density of (?) (color online) 





2DLi2N 


2DLi2N-fH 


2DLi2N+2H 


2DLi2N-fLi 


2DLi2N+2Li 


Eft (eV) 


10.36 


15.98 


16.88 


13.04 


15.37 


Eadd (eV) 




5.80 


0.90 


2.86 


2.29 


djv-_tf (A) 




1.046 


1.157 






djv-ii (A) 








1.88 


1.99 


AZN-Li (A) 


0.0 


0.129 


0.0 


0.187 


0.0 


TABLE I. Gomparison 


of binding and addition energies and 




structure parameters of two-dimensional plane 2DLi2N. 






2DLi2N 


2DLi2N-fH 


2DLi2N-f2H 


2DLi2N-(-Li 


2DLi2N-f2Li 


Et (eV) 


0.54 


0.29 


1.08 


0.21 


0.01 


Eadd (eV) 




5.55 


1.44 


2.11 


1.74 


dN~H (A) 




1.05 


1.21 






djv-ii (A) 








1.91 


2.00 


AzN-Li (A) 


0.49 


0.60 


0.84 


0.30 


0.22 



TABLE II. Gomparison of total energy (with reference to 
plane structure) and structure parameters of two-dimensional 
low-buckled 2DLi2N. 




FIG. 3. Bandstructures and partial density of states for PL 
2DLi2N, 2DLi2N+H and 2DLi2N+2H. Details in text. 
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FIG. 4. Bandstructures and partial density of states for PL 
2DLi2N, 2DLi2N+Li and 2DLi2N+2Li. Details in text. 
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FIG. 5. Bandstructures and partial density of states for LB 
2DLi2N, 2DLi2N+H and 2DLi2N+2H. Details in text. 
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FIG. 6. Bandstructures and partial density of states for LB 
2DLi2N, 2DLi2N+Li and 2DLi2N+2Li. Details in text. 
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